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Abstract. 

Based on a combination of cluster dynamical mean field theory (DMFT) and 
density functional calculations, we calculated the angle- integrated spectral density in 
the layered s = 1/2 quantum magnet TiOCl. The agreement with recent photoemission 
and oxygen K-edge X-ray absorption spectroscopy experiments is found to be good. 
The improvement achieved with this calculation with respect to previous single-site 
DMFT calculations is an indication of the correlated nature and low-dimensionality of 
TiOCl. 
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1. Introduction 

The low-dimensional quantum spin system TiOCl has received a lot of attention in 
recent years due to its anomalous behavior in a wide range of temperatures. It consists 
of Ti-0 bilayers in the afe-plane separated by layers of Cl~ ions stacked along the 
crystallographic c-axis (see FigED). Ab initio density functional calculations^ El E] 
showed that the system at room temperature can be described in terms of spin- 1/2 
(Ti 3+ ) Heisenberg chains running along the crystallographic 6-axis, alongwith small but 
non-negligible inter-chain couplings. Upon varying temperature, TiOCl undergoes two 
successive phase transitions, one of second order nature at T c2 =91K and one of first order 
nature at T c i=67K to a spin-Peierls dimerized state [H 0]. The nature of the transition 
at T c2 =91K has been discussed within various scenarios of orbital fluctuations^ [5], [6], 
frustration of interchain interactions^ [8], the role of phononsfH [HH [TT], spin and 
correlation effect s[T2l [13] and is still debatable. Even the behavior of TiOCl in the 
high-temperature phase has not been understood in a satisfactory way [14]. 

TiOCl can be viewed as a Mott insulator where the insulator gap is driven by 
electron correlation. Attempts to describe the influence of correlation on the properties 
of TiOCl have been performed by considering the local density approximation (LDA)+U 
approach[U [2J, [3] as well as the more elaborate single-site LDA+dynamical mean 
field theory (DMFT) approach[12], [13] where dynamical fluctuations - absent within 
the LDA+U approach - are considered. Nevertheless, comparison of the spectral 
function obtained from single-site LDA+DMFT calculations with photoemission (PES) 
experiments show only a moderate agreement [H]. While the width of the computed 
spectral function is in accordance with the PES results, the magnitude of the optical 
correlation gap was highly underestimated. Also the agreement of the line shape is not 
very satisfactory. 

The effective low-dimensionality of the material suggests that correlation-driven 
inter-site fluctuations may be important for the proper description of TiOCl. In fact, 
it was shown in a recent work[15j on a one-dimensional extended Hubbard model that 
consideration of two-site clusters within the cluster extension of DMFT (cluster-DMFT) 
leads to a satisfactory description of the charge gap. One may therefore expect for TiOCl 
-which is effectively a low-dimensional system- that an improvement upon previous 
single-site LDA+DMFT results may be achieved by including nonlocal correlation 
effects. Demonstration of such a property in a real material of the complexity as TiOCl 
has not been done previously. In order to analyze this proposal, we perform in the 
following calculations within cluster-DMFT where Ti-Ti pairs/dimers are taken as the 
basic units instead of a single site. 

We use the LDA Wannier functions obtained by an N-th order muffin-tin orbital- 
based (NMTO) downfolding method[16] to construct a Hubbard Hamiltonian which we 
solve within the cluster-DMFT approximation. 
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2. Crystal Structure and Energy Levels 

The octahedral environment of Ti [TiC^C^] in TiOCl splits the five degenerate d levels 
into t 2g and e g blocks. Since the Ti 3+ ion is in a 3d 1 configuration, the t 2g states are 
l/6th filled[U[2] and the TiOCl system can be described by a low-energy multiband, t 2g 
Hubbard Hamiltonian [12] . The [Ti04Cl2] octahedra are quite distorted, the t 2g states, 
therefore, show further splittings into lower d xy and higher gL = ^(d xz — d yz ) and d + 
= -^{d xz + d yz ) orbitals, with the z axis pointing along the crystallographic a-axis and 
x and y axes rotated 45° with respect to the crystallographic 6-axis and c-axis. With 
this choice of axes one has the usual convention of d x 2_ y 2, and d z i pointing towards O 
and CI neighbors. This however leads to a non-diagonal form of the Hamiltonian. The 
LDA density matrix (M) for TiOCl has then matrix elements between the d xz and d yz 
orbitals which are nonzero < d xz \M\d yz >= < d yz \M\d xz > ^ 0. A representation of 
the density of states in this basis reflects only the information of the diagonal matrix 
elements < d xz \M\d xz > = < d yz \M\d yz > (see Fig. 2 of Ref. [2]). Diagonalization of 
this matrix provides the corresponding representation into gL and d + . If one chooses 
instead a coordinate system with z = a, x = b, y = c then the Hamiltonian is diagonal 
with d x 2_ y 2, d yz and d xz forming the t 2g block and d xy , d z 2 forming the e g block[17j. 

In the following we will work in the coordinate system with d xy , cL and d + as t 2g 
states. The t 2g crystal field splittings calculated within NMTO are Ai(d xy <-> gL) = 
0.29 eV and A 2 (d xy <-> d + ) = 0.56 eV which are in reasonable agreement with infrared 
absorption spectroscopy measurements [18] where A 2 = 0.65 eV and also comparable to 
recent cluster calculations for the TiO/tCkj octahedronpTS] Ai = 0.25eV and A 2 = 0.69 
eV. 



3. LDA+DMFT calculations 

Following our previous work in Ref. [12], we consider as the starting point for 
the description of TiOCl, the LDA-NMTO Hamiltonian, H m ^, to which direct and 
exchange terms of the screened onsite Coulomb interaction U mm ' and J mm ' of Hubbard- 
Hund type are added [T9]. 

H = '^liR m ,R'm l ,(T ^Rm%'m' ^RmpR'm' + 2 51 ^rnm'^R rn P J Rm' W 

R,m,m',cr 

+ 2 S (Umm' — Jmm') ^Rm^Rm' 
R,m^m' ,cr 

Here, c R ^ m (c Rm ) denotes the creation (annihilation) operator for an electron at site R 
in orbital m with spin a, and n Rm = c Rm c Rm is the particle number operator. U mm < 
and Jmm' are parametrized as follows: U mm = U is the Coulomb repulsion between 
electrons in the same orbital, U m ^ m i = U — 2 J is the average repulsion, and J m ^m' = J 
is the Hund's rule coupling. We assume the double counting correction to be orbital 
independent within the t 2g states, thus resulting in a simple shift of the chemical 
potential. The choice of the local orbitals ( Ti-t 2g Wannier functions) is done via 
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Figure 1. Left Panel: Crystal structure of TiOCl. Big spheres represent Ti atoms and 
small spheres represent O and CI atoms. The smallest box shows the supercell used for 
the LDA calculation. Two nearest-neighbor dark and light shaded Ti atoms positioned 
along the crystallographic 6-axis form the pairs. Right Panel: Network of Ti atoms 
within the supercell projected in ab plane. Two different sizes of atoms correspond 
to Ti atoms belonging to two different layers of the bilayer. Thick and thin bonds 
represent the intra and inter-dimer bonds respectively in the supercell. 

the NMTO-downfolding technique [16], where all the partial waves other than Ti— d xy , 
Ti— d xz and Ti— d yz are downfolded (these notations refer to the choice of coordinate 
system as explained previously). As shown in Ref. [12] and earlier communications 
[20], the recent implementation of LDA+DMFT[2T] allows to solve the many-body 
Hamiltonian, including all off-diagonal elements in the orbital space of the local self- 
energy, S mm / . This has been crucial for TiOCl since the choice of the coordinate system 
as discussed above does not result into a diagonal form of the onsite matrix. 

Our starting point is the high-temperature crystal structure. Since the low- 
temperature crystal structure of TiOCl shows a doubling of the cell along the 
crystallographic b natural choice of Ti pairs are those along the b axis. For 
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our cluster calculation we have therefore carried out a supercell calculation with the 
unit cell doubled along the b— axis. This results into four Ti atoms in the unit cell, 
marked as 1, 2, 3 and 4 in Fig. [U with two Ti-Ti pairs located in the upper and lower 
layer of a given bilayer. Our self-energy will have accordingly off-diagonal elements in 
the orbital space as well as in the site space. The subscript M in the self-energy Emm' is 
defined as a composite site-orbital index. Such formulation has been already successfully 
applied in the case of VO2 and Ti 2 03|22, 23j. 

For the present problem, a 12 x 12 block-diagonal self-energy matrix, Emm' is 
therefore constructed where M is the composite index (m, i c ) with m denoting the orbital 
index, 1,2,3 for a t 2g basis and i c = (1,2) and (3,4) denote the intradimer site indices for 
two Ti-Ti pairs in the unit cell. We neglect the interdimer correlation connecting the two 
pairs (1,2) and (3,4), and set the interdimer components of the self-energy in the present 
calculation to zero resulting in a block diagonal form of the self-energy. This may be a 
reasonable approximation considering the fact that the effective intradimer interaction 
between 1 and 2, or 3 and 4 is about an order of magnitude larger compared to interdimer 
interactions connecting two pairs belonging to different layers in this system[2]. With 
this approximation, the self-energy takes the form: 
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where E 11 and E 22 denote the on-site self-energy corresponding to sites 1 and 2 
within a pair. Each of these matrices, is therefore a 3 x 3 matrix. E 12 (E 21 ) gives 
the intersite, intradimer component of the self-energy. Note that the presence of the 
intersite component of the self-energy, E 12 (S 21 ) gives rise to some[24] ^-dependence of 
E, as expected for a cluster-DMFT calculation. 

We further assume that the dimers belonging to two different layers of the bilayer 
are similar, therefore the upper and lower sub-blocks of the matrix given by (1) are 
identical i.e. E 33 = E 11 , E 44 = E 22 and so on. It is therefore enough to work with a 6 x 
6 block of this self-energy in the DMFT self-consistency condition and then construct 
the full 12 x 12 block from the 6x6 block to which H^ffi is added to generate the 
Green's function. One may note that, with such a choice, the bond ordered ground 
state is allowed in the calculations, although one needs to include the electron-phonon 
interaction in the Hamiltonian in order to drive that state. For a single chain of Ti 
atoms running along the b-axis, the dimers may be formed in even or odd bonds. 
For a bilayer, as is the present case, there are however four possible arrangements |8j, 
with even/even, odd/odd, odd/even and even/odd bonds in upper/lower layers. With 
our present choice of supercell it is not possible to distinguish between these cases. 
Consideration of such scenarios would require to have a supercell which is four times 
enlarged along the 6-axis. A second possible dimer pattern scenario, which is also not 
included in the present calculation, is the in-phase/out-of-phase dimer arrangements of 
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neighboring chains within a layer. Such arrangements can be distinguished in a supercell 
calculation, which in addition to doubling along 6-axis is twice enlarged along the a-axis. 
Both calculations are presently computationally too expensive. 

The 6x6 impurity problem is solved by a numerically exact Quantum Monte 
Carlo (QMC) scheme. Within the given computational resources, we could reach a 
temperature of 1400 K with 10 6 QMC sweeps and 68 time slices. U and J values were 
chosen to be 4 eV and 0.7 eV respectively. Whereas J = 0.7 eV is a generally accepted 
value for an early transition- metal oxide (hardly changed from its atomic value), the 
precise determination of U is a delicate issue. Ab-initio techniques such as constrained 
LDA give only rough estimates. For an early transition metal like Ti, U is expected to 
be between 3 to 5 eV[2]. We carried in the past calculations [32] with different choices 
of U in terms of a single-site DMFT calculation and fixed U to 4.0 eV since it provided 
the best possible spectra in terms of band gap. We consider here the same U and 
J values as above so as to compare the single-site and cluster-DMFT results. The 
maximum entropy method [25] has been employed for the analytic continuation of the 
Green's function from the imaginary to the real axis for the calculation of the spectral 
functions. 

Experimentally, the electron removal and addition spectra were probed by angle- 
integrated photoemission spectroscopy (PES) and x-ray absorption spectroscopy (XAS), 
respectively. Details of the PES experiment are described in Ref. [T4] . XAS was 
measured at the PM-3 beamline of BESSY (Berlin, Germany) using its MUSTANG 
endstation. Due to finite Ti 3d-0 2p hybridization the oxygen K-edge spectrum can be 
taken as an approximate measure of the Ti 3d electron addition spectrum. The energy 
resolution for PES and XAS amounts to 100 meV and 50 meV, respectively. 

4. Results 

Fig. [2] shows the Ti t 2 g-dominated total spectral function computed with the 
NMTO+cluster-DMFT method described above (full line) in comparison with the angle 
integrated photoemission data (blue solid dots) and oxygen K-edge XAS measurements 
(red open dots) in the range of energies between -5eV and 7eV. Note that the 
experimental XAS spectrum shows both the t 2g and e g manifolds. The insulating 
behavior of TiOCl is correctly described by this calculation with a charge gap of about 
1.1 eV, in rough agreement with the observed 2 eV optical gap[8l 126]. The important 
point to note is the enhancement of the gap value compared to single-site DMFT result 
which was about 0.3 eV[12]. This inevitably points towards the importance of the inter- 
site fluctuations for a good description of a low-dimensional system like TiOCl. 

The calculated spectral weight distribution below and above the chemical potential 
jj, i.e. the lower and upper Hubbard bands compare reasonably well with the PES and 
oxygen K-edge absorption data, respectively, though one should keep in mind that the 
temperature used in the theoretical calculation was rather high compared to that of the 
experiments. We observe that the spectral function below [i is dominated totally by 
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Figure 2. Comparison of the cluster-DMFT spectral function (full line) with 
angle integrated photoemission data (blue solid dots) and oxygen K-edge absorption 
spectrum (red open dots). The experimental spectra have been aligned in energy 
position and intensity to the theoretical spectrum. 

the (i^-like contribution with a 99% occupancy, consistent with polarization-dependent 
optical spectroscopy [8] and ARPES[H] experiments. Note that our calculations do not 
show any shape for the observed empty e g states since these bands were not explicitly 
considered in the cluster-DMFT calculations. They were downfolded and included only 
as tails of the t 2g NMTO Wannier orbitals and may be responsible for the slightly larger 
width of the calculated upper Hubbard band compared to the t<i g peak of the XAS 
spectrum. 

In Fig. [U we present a comparison of the cluster-DMFT with the single-site 
DMFT resultspE2]. Both calculations were performed within the same NMTO basis 
and considering the same QMC impurity solver. We notice improvement of the present 
results with respect to the single-site DMFT results This may emphasize once more the 
importance of including dynamical Ti-Ti intersite correlations for the description of the 
spectral properties of TiOCl. For comparison, we also show the spectra obtained from 
LDA+U calculations which has a much narrower width as has been already stressed in 
Ref. [H]. 

While the accordance between theory and experiment turn out to be reasonably 
good as we have shown above, there are a few sources of improvement, which are all 
related to the computer feasibility of the calculations and are out of the scope of the 
present work, i) The presented results were performed at quite high temperatures 
and with moderate numbers of QMC sweeps and time-slices steps regarding the 
QMC impurity solver. Calculations at lower temperatures are expected to provide 
a better resolution of the results. In particular, it will be necessary to check the 
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Figure 3. Comparison of cluster-DMFT, single-site DMFT spectral functions and 
LDA+U resultspj] with photoemission. To facilitate better comparison between 
experimental and theoretical lineshapes an inelastic background has been subtracted 
from the photoemission spectrum. All spectra are aligned to the same first order 
moment and scaled to the same integrated weight. The absolute energy scale refers to 
that of experiment. 



temperature dependence of the spectra which is at the moment beyond our capability, 
ii) Consideration of larger supercells is expected to distinguish between the odd/even 
and even/even arrangements as well as in-phase/out-of-phase patterns and may provide 
a clue of the possible bilayer frustration present in this system, iii) In order to be both 
computer realistic and manageable, we have used the NMTO ti g Wannier functions for 
a low-energy subset of LDA bands, while in the actual case there are also O and CI— p 
dominated bands and also higher lying Ti-e g states. The influence of the CI and O-p-like 
bands (and also Ti-e 9 bands) has been considered implicitly in the construction of the 
t 2g Wannier functions which take care of the hybridization effects coming from CI and 
O-p-like bands (see Fig. 3 in Ref. [2]), but they have not been considered explicitly. 
In principle, one should therefore carry out calculations involving O and Cl-p bands 
too, which is an extremely computer intensive job and not possible to carry out in the 
present framework. 

5. Conclusions 

- By means of cluster-DMFT calculations implemented in the NMTO Wannier function 
basis, we have shown that the Ti-Ti intersite correlations in TiOCl play an important 
role for the proper description of photoemission and O K-edge XAS spectra. Although 
our calculations were limited by the computer feasibility, the improvement of the value 
of the optical gap and the overall comparison of the theoretical spectra with the 
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experimental one indicates the significant role of the off-site, intra-dimer correlations. 
The present results suggest that the nature of fluctuations observed in a large variety 
of experiments [5l [TO] may be governed by correlated Ti-Ti dimers, existing up to rather 
high temperatures. 
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